Negative-coupling resonances in 
pump-coupled lasers 



T.W. Carr a '*, M.L. Taylor a 

^Department of Mathematics 
Southern Methodist University 
Dallas, TX 75275-0156 
tel: 214-768-3460 fax: 214-768-2355 

LB. Schwartz h,x 

b Nonlinear Dynamical Systems Section, Code 6792 
Plasma Physics Division 
Naval Research Laboratory 
Washington, DC 20375 
tel: 202-404-8359 fax: 202-767-0631 



Abstract 

We consider coupled lasers, where the intensity deviations from the steady state, 
modulate the pump of the other lasers. Most of our results are for two lasers where 
the coupling constants are of opposite sign. This leads to a Hopf bifurcation to 
periodic output for weak coupling. As the magnitude of the coupling constants is 
increased (negatively) we observe novel amplitude effects such as a weak coupling 
resonance peak and, strong coupling subharmonic resonances and chaos. In the 
weak coupling regime the output is predicted by a set of slow evolution amplitude 
equations. Pulsating solutions in the strong coupling limit are described by discrete 
map derived from the original model. 
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1 Introduction 



In recent work, we presented experimental and simulation results for two cou- 
pled lasers [1] with asymmetric coupling. That is, the coupling strength from 
laser- 1 to laser-2 was kept fixed, while the coupling strength from laser-2 to 
laser- 1 was used as a control parameter. In this paper, we present a more 
theoretical exploration of the dynamics that result from this coupling config- 
uration. Each laser is tuned such that it emits a stable constant light output. 
Light-intensity deviations from the steady state are converted to an electronic 
signal that controls the pump strength of the other laser. Our work in [1] 
considered asymmetric coupling and, more specifically, the effect of delaying 
the coupling signal from one laser to another. The present paper is our first 
theoretical analysis of two pump-coupled lasers with asymmetric coupling, but 
without delay (analysis of the case with delay will be presented in a future 
manuscript). However, we do invert one of the electronic coupling signals such 
that the effective coupling constant is negative; for harmonic signals with de- 
lay coupling by half the period, both lead to the same phase shift. Thus, the 
present paper serves as a prelude to a future study of two pump-coupled lasers 
with same-sign delay coupling. 

For very weak coupling, both lasers remain at steady state. As the coupling 
is increased, but still small, there is a Hopf bifurcation to oscillatory output. 
In the weak-coupling regime we also observe and describe a resonance peak 
where the amplitude of both lasers becomes large over a small interval of the 
coupling parameter; to our knowledge, this phenomenon has not previously 
been reported. For strong coupling, the oscillations of one laser remain small 
and nearly harmonic while the other laser exhibits pulsating output. Period- 
doubling bifurcations to chaos and complex subharmonic resonances also exist 
throughout the parameter regime. We combine both weakly- and strongly- 
nonlinear asymptotic methods to describe the output in the case of strong 
coupling. 

We consider two class-B lasers [2,3] modeled by rate equations as 



^± = ( Dj -l)I j: j = 1,2 (1) 

^=^-{1 + 1^}, (2) 

where Ij is intensity and Dj is the population inversion of each laser. Dimen- 
sionless time t is measured with respect to the cavity-decay time k , t = k t r , 
where t r is real time. The parameters are 

«»_£, A=^f P , (3) 
ko Icko 
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where e 2 is a ratio of the inversion-decay time, j c , to the cavity-decay time, ko, 
and A is proportional to the pump P (for notational clarity we have suppressed 
the subscript j on the parameters in the definitions of e and A). To facilitate 
further analysis, we define new variables for the deviations from the non-zero 
steady state (CW output) [4] Dj = 1, Ij = Aj — 1 as 



Ij = I j0 (l +Vj), Dj = 1 + tj^IjoXj, S = €iy/I 10 t. (4) 

Our goal is to investigate the effects of coupling through the pump with 
Aj = A j0 + I j0 5 k yk- (5) 

Thus, we feed the intensity deviations y k = (Ik — ho)/ho from the CW output 
of laser k to the pump of laser j; the strength of the coupling is controlled 
by 5k- The pump coupling scheme allows for easy electronic control of the 
feedback signal. 

Finally, we assume that the decay constants of the two lasers are related by 
e 2 = ei^=(l + €ia). The new rescaled equations are 



— = x 1 (l + y 1 ), 

dx\ , . . . 

= -Vi - exi(ai + by!) + d 2 y 2 , 

^ = px 2 (l+y 2 ), 

^- = P[-y2-ePx 2 (a 2 + by 2 )+5 1 y 1 ], (6) 



where 



^ + \/7lo(i + ho) , rr~ Aft i i C7\ 

ai = — — , a 2 = , b = y ho and p = 1 + ea. (7) 

V ho ho 



For notational convenience we have let s — + t and dropped the subscript on 
ei(ei — > e). We mention that Eq. (6) is similar to the coupled laser equations 
studied by Erneux and Mandel [5] to investigate antiphase (splayphase) dy- 
namics in lasers. However, antiphase dynamics require global coupling that 
would correspond to the symmetric case of 5 2 = Si in our model. 

A main point of interest in the study of coupled oscillators in general is their 
degree of synchronization. This implies a focus on the phase- and frequency- 
locking characteristics of the oscillators. Thus, many investigations focus on 
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coupled-phase oscillators (see [6] and [7] for reviews and extensive bibliogra- 
phies). Consideration of just the phase relationships between the oscillators 
is often based upon considering limit-cycle oscillators with weak coupling. In 
that case, each oscillator's amplitude is fixed to that of the limit cycle and 
only the phase remains a dynamical variable. However, limit-cycle oscillators 
with strong coupling can exhibit amplitude instabilities leading to amplitude 
death and other novel phenomena [6]. 

Class-B lasers, which include such common lasers as semiconductor, YAG, and 
CO2 lasers, are not limit cycle oscillators, but, rather, are perturbed conserva- 
tive systems [8]. (The underlying form of the perturbed conservative system, 
Eqs. (6) with e = 0, has also been used in population dynamics models [9].) 
Thus, the amplitude is not fixed by a limit cycle and remains an important 
dynamical variable. This has been demonstrated in laser systems coupled by 
mutual injection [10] and overlapping evanescent fields [11], or by multimode 
lasers with coupled modes [12,19] to name just a few. Under certain condi- 
tions, phase-only equations can be derived that describe the behavior of the 
coupled laser systems [11,20,21]. However, in general, the amplitude cannot be 
adiabatically removed and amplitude instabilities can dominate the observed 
dynamics. In is interesting to note, however, that a time-dependent phase is 
sometimes the drive leading for the laser's observed amplitude instability [22]. 

The coupled laser equations, Eqs. (6), have all real coefficients. If the lasers 
were coupled directly through their electric fields (referred to as "coherent 
coupling"), such as in evanescent or injection coupling, then there would be a 
complex detuning parameter or coupling coefficient. In Eqs. (6) the lasers are 
coupled through their real intensities (referred to as "incoherent coupling") 
such that the differences in the laser's optical frequency do not affect the 
systems dynamics. 

In the next section, we give an overview of the laser system's behavior as 
the coupling is increased. We begin with the linear-stability analysis of the 
CW steady state (x, y) = (0, 0) and find that there are two possible Hopf 
bifurcations to oscillatory output, one for large 0(1) coupling, and one for 
small 0(e) coupling; we focus the rest of our analysis on the latter and continue 
our overview by presenting results of numerical simulations over the full range 
of coupling strengths. In Sec. 3 we analyze the oscillatory solutions for weak 
coupling. We also show how the results in this parameter regime extend to the 
case of three or more lasers. In Sec. 4 we consider large coupling and combine 
the method of multiple scales and matched asymptotics [13] to derive a map 
that describes the coexisting small- and large-amplitude solutions. Finally, in 
Sec. 5 we discuss and summarize our results. 
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2 Bifurcations for negative coupling 



In the new variables, the CW state is given by x = y = 0. The linear stability 
of the CW state is governed by the characteristic equation 

[A(A + eai) + 1] [A (A + ea 2 (3 2 ) + (3 2 } - (3 2 SiS 2 = 0. (8) 



If both 5j = 0, as expected we find that each laser is a damped oscillator. For 
5j ^ we study Eq. (8) for small e € 1. Keeping Si as a fixed parameter 
and varying S 2 we find that there are two Hopf bifurcations. If S 2 < and \S 2 \ 
increases, then the condition for a Hopf bifurcation is 

SJ 2 + e 2 [ ai a 2 + Aa 2 . — ._ ] + 0(e 3 ) = 0. (9) 
(oi + a 2 ) 1 

If 5 2 > and \5 2 \ increases the Hopf condition is 

5 1 5 2 = l + 0(e 2 ). (10) 



The second condition, Eq. (10), indicates that a Hopf bifurcation occurs when 
there is strong coupling between the lasers, 5i5 2 = 0(1). We are interested in 
the Hopf bifurcation that occurs for weak coupling that is described by the 
first condition, Eq. (9) (this is the relevent case when the problem is extended 
to include delayed coupling). In this case, 5 2 h = 0(-e 2 /5 1 ) < 0, that is, the 
coupling from laser-2 to laser- 1 is negative. 

In Fig. A.l we show the amplitude of the periodic solutions that emerge from 
the Hopf bifurcation point of Eq. (9). As the magnitude of the coupling con- 
stant (\S 2 \, S 2 < 0) is increased, the Hopf bifurcation leads to small-amplitude 
periodic solutions. However, for small coupling there is a strong resonance ef- 
fect where the amplitudes become 0(1). In Fig. A.l this appears as a narrow 
spike in the amplitude. We show a close-up of the amplitude resonance in 
Fig. A. 4 (calculated at different parameter values). Both before and after, 
the amplitude is small and nearly harmonic, as would be expected for weak 
coupling. However, during the resonance the amplitude is pulsing. 

As \5 2 1 is increased, the coupled system behaves similar to a periodically mod- 
ulated laser [4,14]. The intensity of laser- 1 increases and becomes pulsating 
(see Fig. A. 2a) because the effective modulation signal from laser-2 becomes 
stronger. On the other hand, because Si is fixed and small, laser-2 receives 
only a weak signal from laser- 1 and remains nearly harmonic (see Fig. A. 2b). 
For larger coupling, the periodic solutions exhibit a period-doubling sequence 
to chaos; the inversion for both laser- 1 and -2 after the first period-doubling 
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bifurcation is shown in Fig. A. 3b. We mention also that for different parameter 
values the original branch of periodic solutions may remain completely stable 
and not exhibit further bifurcations. 

Coexisting with the primary branch of periodic solutions are subharmonic 
resonances that appear through saddle-node bifurcations. These also exhibit 
period-doubling bifurcations for increasing coupling. In Fig. A. 3c we see that 
just after the primary saddle- node bifurcation the periods of the oscillations 
are in a 2:3 ratio, with 2 maximum of laser- 1 for every 3 of laser-2. 



3 Weak-coupling resonance 



3.1 Two lasers 



We now describe the periodic solutions that emerge from the Hopf bifur- 
cation located by Eq. (9). We use the standard method of multiple time- 
scales [13] approach and thus only summarize the results. From the linear- 
stability analysis, we know that solutions decay on an 0(e) time scale. This 
suggests that we introduce the slow time T = et, such that x = x(t,T) 
(similarly for y) and time derivatives become ^ = + l 2 -§f- We analyze 
the nonlinear problem using perturbation expansions in powers of e 1 / 2 , e.g., 
Xj(t) = e l l 2 Xji(t,T) + eXj 2 (t,T) + . . ., Finally, we assume that the coupling 
constants are small and let 5j = edj. 

At the leading order, 0(e 1//2 ), we obtain the solutions 

y n (t, T) = A j (T)e it + c.c, x(t, T) = iA j (T)e it + c.c, (11) 



which exhibit oscillations with radial frequency 1 on the t time scale. To find 
the slow evolution of Aj(T) we must continue the analysis to 0(e 3 / 2 ). Then, to 
prevent the appearance of unbounded secular terms, we determine "solvability 
conditions" for the Aj (T) that are given by 



dAi 1 1 .. , l9 , 1 . , , 

-j-L = --a 1 A 1 - -i\A 1 \ 2 A 1 - -zd 2 A 2 , (12) 

dAo 1 . 1 , . ,9 . 1 , . . /N 

— - = --a 2 A 2 - -i\A 2 \ 2 A 2 - -id x A x + iaA 2 . 13 
al I o I 

To analyze these equation we let Aj(T) = Rj{T)e %63 ^ and consider the phase 
difference ip = 9 2 — Q\ to obtain 
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^ = -^a 1 R 1 + ^sm(^)d 2 R 2 , (14) 
( ^- = -^a 2 R 2 -^sm(^)d 1 R 1 , (15) 
^ = -l(Rl -R\)-- cos iP(d^ - dA + a. (16) 

The leading order, solutions are t = 27r periodic if the amplitudes and phase 
are constant with respect to the T time scale (derivatives with respect to T 
are zero). This determines the bifurcation equation for the amplitudes Rj and 
the phase difference ip as 



t = -9^, where A 1 = (l + — )(ai + a 2 ) 2 , and A 2 = 1 + 

A 2 a x a 2 ' a\d\ 



and 



#i = 1-^ 2 , cos ^ = (1 + ——), (18) 

aicsi d\d 2 



where we have set a = to simplify the discussion. For Ri to be positive 
in Eq. (18), d\ and d 2 must have opposite signs, while the Hopf bifurcation 
point is determined by taking R 2 — * in Eq. (17) to obtain Ai = 0; both 
conditions are consistent with the linear stability results in Eq. (9). We define 
the value at which the Hopf bifurcation occurs to be S 2 h = zd 2 H- For d 2 > d 2 u, 
Eq. (18) describes a supercritical bifurcation to stable periodic solutions; this is 
consistent with the numerical bifurcation diagram in Fig. A.l. Finally, because 
dip/dT = 0, the laser oscillations are phase locked with the phase difference 
described by Eq. (18), and the frequency for x and y is 

2 V \|A 2 | a 1 + a 2 J 

An important result of this paper comes from an examination of A2 in Eq. 
(17). Specifically, the bifurcation equation is singular when A 2 = or 

d 2 = d 2S = -—d 1 . (20) 
a 2 



If d 2 ,s < d 2 H, then the singularity occurs before the Hopf bifurcation when 
the CW steady-state is still stable. Thus, in this case, the singularity is not 
seen and does not affect the amplitude of the bifurcating periodic solutions. 
However, if d 2 $ > d 2 n, then near the bifurcation point the amplitude of the 
oscillations becomes very large corresponding to a resonance. The resonance 
can be understood as a balance between an effective negative damping due to 
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the coupling term, and the self damping. That is, the ratio d 2 /a\, which is the 
relative negative camping to the self damping in laser- 1, is equivalent to d\ja1 
(modulus the negative sign), the relative negative damping to self damping in 
laser-2. The net result is that the coupling terms provide an effective negative- 
damping that cancels with the lasers self-damping and, hence, a resonance 
effect. 

The negative-coupling resonance when d 2 = d 2 s is demonstrated in Fig. A. 4a. 
The solid line is the result of our analytical bifurcation curve given by Eq. (17), 
while the + are data from numerical simulation; the analytical and numerical 
results are in excellent agreement. In the vicinity of 5 2 = S 2 s the amplitude of 
the periodic oscillations become 0(1), whereas we would normally expect the 
amplitude to remain 0{e 1 ^ 2 ). 

Comparing Fig. A. 4a and Fig. A. 5, we see that the maximum amplitude, 
when 5 2 = 5 2 $, depends on the parameters. However, the bifurcation equation 
is singular at 5 2 s and does not give a value for the maximum. The bifurca- 
tion equation can be improved by tuning the resonance closer to the Hopf 
bifurcation point with Si = ea 2 + 0(e'l 2 ) and 5 2 = —ea>i + 0(e 3 / 2 ) and contin- 
uing the perturbation analysis to 0(e 2 ). Unfortunately, the analysis become 
algebraically difficult and we have not pushed through to its conclusion. 

During the resonance both lasers become pulsating. Pulsating solutions are 
not well described by the weakly-nonlinear analysis of the present section. In 
Appendix A we consider pulsating lasers and again locate the resonance peak 
at S 2 s- We discuss this further in the paper's final discussion section. 

3.2 Three (or more) lasers 

The resonance spike can also be found in three or more lasers. In general, the 
amplitudes of the periodic solutions near the Hopf bifurcation are described 
by coupled Stuart-Landau equations of the form 



As written, the coupling coefficients are completely general and could be cho- 
sen to give global coupling, djk = d, nearest-neighbor coupling, dj t k ^ for 
only k — j + 1, k — j — 1, or any other coupling configuration. Coupled 
algebraic equations for the amplitudes Rj can then be found with the sub- 
stitution Aj = Rj exp(iOj). An amplitude resonance occurs when there is a 
vanishing denominator in the equation for any one of the isolated amplitudes 
Rj = g(Rk), k 7^ j. As with two lasers, one of the coupling constants must be 



dAj 




i\Aj\ 2 Aj--i ]T d jk A k , j = l...N. (21) 



dT 



k=i,k^j 
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negative to produce the resonance. However, obtaining an explicit solution for 
one of the laser amplitudes, even in the case of only three lasers, is extremely 
difficult in all but the most trivial cases. 

In contrast, demonstrating the resonance effect numerically requires only some 
experimentation and we show one result in Figs. A. 6 and A. 7. In Fig. A. 6 
we show the amplitude of each laser as a function of one of the coupling 
parameters. Specifically, we fix the coupling of laser-3 into laser- 1 and laser-2 
as dis = (I23 = 1.3 and the coupling of laser- 1 into laser-2 and laser-3 as d 2 \ = 
d 3 i = 3. The coupling of laser-2 into laser-3 is positive with size d 32 = \d 2 \, 
while the coupling of laser-2 into laser-1 is negative with d 12 — d 2 < 0. We use 
d 2 as the control parameter. 

As \d 2 \ is increased, both laser-1 and laser-2 show two resonance peaks, while 
for laser-3 there is only one. However, in Fig. A. 7 we see that the branch of 
solutions is not monotonic in \d 2 \. As the branch of solutions is followed from 
the Hopf bifurcation point, the resonance peak for larger \d 2 \ ~ 3 occurs first, 
but only for laser-1 and laser-2. As the branch is followed further, it turns 
at the saddle- node bifurcation (right most in figure). As \d 2 \ decreases, all 
three lasers exhibit a resonance when \d 2 \ ~ 1.75. The branch turns again 
at a saddle-node bifurcation (left most) to then increase without any further 
resonances. 

As it happens, the periodic solutions are unstable on the branch of solutions 
with the lower resonance peak \d 2 \ ~ 1.75. Thus, between the two saddle-node 
bifurcations there are two stable solutions: the primary branch originating 
from the Hopf bifurcation that exhibits a resonance for laser-1 and laser-2 
and terminates at the larger (right) saddle-node bifurcation, and the small- 
amplitude branch that appears at the lower (left) saddle-node bifurcation and 
continues for \d 2 \ > 5. Thus, in the vicinity of the resonance when \d 2 \ ~ 3 
and both stable solutions coexist, initial conditions will determine whether 
the large-amplitude resonant solutions or the small-amplitude solutions are 
exhibited 

Finally, the period of the oscillations shows a sharp peak at each of the am- 
plitude resonances. The result is analogous to that of the peak in the period 
for two lasers as shown in Fig. A. 4b. 



4 Strong coupling 

For "large" values of the coupling, when 5 2 = 0(1), the intensity of laser-1 
becomes pulsating, while the oscillations of laser-2 remain small and nearly 
harmonic (see Fig. A. 2). We derive an iterated map to describe the oscillations 
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when 82 = 0(1). Fixed points of the resulting map correspond to periodic 
solutions of Eq. (6). Our results are summarized in Fig. A. 8, where we compare 
the amplitudes and period to those obtained from numerical simulation. 

To construct the map we take advantage of the fact that the intensity of 
laser- 1 has two distinct regimes: during the pulse when yi >> 1, and a long 
interval of time when yi pa — 1. For a single pulsating laser, the solutions to 
Eq. (6) have described using an iterated map constructed with the method 
of matched asymptotics [4,14]; we will use the same approach here and so 
will just summarize our results. We will first find the "outer" solutions to 
Eqs. (6) with the approximation y ± pa — 1 ( see Fig. A. 2a from t to ti). We 
will then reanalyze the coupled system with an "inner" or "boundary- layer" 
approximation y\ >> 1 (see Fig. A. 2a from t\ to t 2 ). The typical next step 
is to match the inner and outer solutions to form a composite solution over 
the whole period. However, we are interested in the dynamics from one pulse 
to the next. Hence, we will simply patch the solutions together to form an 
iterated map. 

As described above, the pulsations of laser- 1 define the inner and outer regimes. 
However, laser-2 continues to exhibit small-amplitude, nearly-harmonic oscil- 
lations. Thus, in each regime we will use the method of multiple scales to 
describe the oscillations of laser-2. 



4-1 Re-supply of the inversion, y\ pa —1 

We first consider the outer regime when y\ pa — 1. We define t = to as the time 
of completion of a previous pulse, when y x — and the inversion x\ is at its 
minimum (see Fig. A. 2). The end of the outer regime will be defined to be 
when the intensity increases from yi — — 1 back to y\ — and the inversion 
x\ is at its maximum. 

We first consider laser-2. When yi pa — 1 the dynamics of laser-2 can be ap- 
proximated as 



We can solve this system using the method of multiple scales as we did in 
Sec. 3.1 under the assumption that Si is small (Si = 0(e)). We find that to 
leading order laser-2 is a weakly damped, nearly harmonic oscillator described 

as 



d]j2 
dt 
dx 2 

~dt 



f3x 2 (l + y 2 ), 



P[-V2 - e(3x 2 (a 2 + by 2 ) - Si], 



(22) 
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x 2 (t)=e 2 ea2 *, [rr 20 cos(cj0) -y 20 sm(u(f))}, <f> = t - t , 

dx 2 



V2(t) 



dt ' 



^=l-\s 1 -^(xl + yl )e-^, (23) 

where (x 2 (t ), 2/2(^0)) — (^20? 2/20) is the state of laser-2 at the end of the 
previous pulse. (As in Sec. 3.1, x 2 and y 2 are 0(e 1 ^ 2 ) such that the next term 
in the solutions in Eqs. (23) would be 0(e).) 

We now examine laser-1 in more detail. With ?/! » -1 we have 

^ = 1 - c(ai - b) Xl + S 2 y 2 , (24) 
at 

which can be integrated to obtain 

1 1 \ 

xi(t) = (xio - -)e-^ + - + 8 2 e-* / e"y 2 (s)ds, (25) 

7 7 J 



to 



where 7 = e(a\ — b). Using the result for y 2 from Eqs. 23 we obtain 



x 1 (f) = (a; 1 o--)e-^ + - 
7 7 

(ay 20 - ujx 2O )(e a0 cos(u;0) - 1) 



cr + u;^ 



+(wy 2 o + ax 2o )e o *sin(w0) , (26) 



where a = (7 — ea 2 /2). We can now use xi to improve our approximation for 
yi by substituting Eq. (26) in the equation for y x in Eq. (6); we then integrate 
to give 

yi (t) = -1 + (1 + y w ) exp ^1 Xl (s)ds^ . (27) 

During the outer regime the inversion grows almost linearly from its minimum 
to maximum values, more precisely, for e 1, x\ ~ (t— 1 ). With the inversion 
re-supplied the laser can then emit a new pulse of light. We define the start 
of the next pulse at t — t\ to be when yi(t\) = 0. Thus, the next pulse begins 
when the integral in the exponential of Eq. (27) is zero, or 

Jx 1 (t)dt = 0. (28) 
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4-2 Pulse regime, y\ >> 1 



The inner regime is defined to be when the intensity is large, y\ >> 1, and 
occurs over a very short interval of time. Specifically if y\ = 0(E), where 
E ^> 1 is related to the energy then X\ = 0(E 1 ^ 2 ) and the width of the 
pulse is 0(1/ E 1 / 2 ) [14]. On the other hand, the oscillations of laser-2 remain 
small, 0(e 1 / 2 ). Thus, we assume that to leading order, laser-2 has no effect on 
laser-1 during the pulse. We can then use the results from [14], in the absence 
of modulation, to describe laser-1. Namely, (i) the end of the pulse, t = t 2 , 
is defined to be when the pulse intensity returns to zero, y(t 2 ) = 0, (ii) the 
width of the pulse is negligible compared to the time in the outer regime, 
(t 2 — t\) « 0, and (iii) the inversion drops from its maximum to minimum 
value with reduction due to damping: 

xi(t 2 ) = -x 1 (t 1 ) + -ebx 1 (t 1 ) 2 . (29) 



(^1(^2) is negative at the minimum so that the additional positive term is a 
reduction in the magnitude of the minimum.) 

The large and narrow (t = 0(1 /E 1 / 2 )) pulse of laser-1 does have a significant 
effect on laser-2. To determine the appropriate inner problem for laser-2, we 
scale the pulse amplitude as y± = 0(E) and stretch time according to t = 
O^/E 1 / 2 ). The coupling is weak with 5\ = 0(e). Finally, we assume that 
E = O^/e 1 ' 2 ) to obtain 



dx 2 

^ = 0. (30) 
dt 

Thus, to leading order y 2 is constant during the pulse while x 2 is given by 
x 2 (t 2 ) - x 2 (h) = 5 1 J Vl (t)dt. (31) 



However, in the pulsing regime we have that dxi/dt ~ —y\ so that 

x 2 (t 2 ) - x 2 (t l ) = -<$! j ^-dt = S^x^h) - Xl (t 2 )]. (32) 
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Thus, for laser-2 we have that 

x 2 (t 2 ) = x 2 (ti) + 5i2x(*i), y(t 2 ) = y(ti) 



(33) 



where we have used Eq. (29) in Eq. (32) and kept only the leading order terms 
(the leading order is 0(e l l 2 ) and we have dropped the 0(e) corrections). The 
net effect is that at the end of the pulse when t = t 2 , the intensity y 2 remains 
unchanged, while x 2 has received a "kick" due to the pulse from laser- 1. 



4-3 Constructing the map 



To construct a map, we "patch" together the results from the outer and inner 
analysis of the previous two sections. 

For laser-2 we initially have (x 2 (to), y 2 (to)) = (x 20 , y 2 o) that in the outer region 
evolves according to Eqs. (23) until t — t±. Then, in the inner region, laser-2 
receives the pulse from laser-1 according to Eq. 33. Thus, we have that 

[(x 2 (t ),y 2 (t )) = (x 20 ,y 20 )] i-> (x 2 {t 1 ),y 2 {t 1 )) ^ (x 2 (t 2 ) , y 2 (t 2 )) . (34) 



The total time from one pulse to the next is t 2 — to- However, because the 
pulse is so short ((^(e 1 / 2 )), we make the approximation that t 2 ~ t\ and define 
the total time as P = t\ — t . Finally, for notation convenience we define the 
intermediate value of the inversion of laser-1 as G(P) = Xi(ti). The map for 
laser-2 is then 



x 2 ^ e -^ ea2P [x 2 cos(cjP) - y 2 sm(uP)] + 5i2G(P), 

y2 e -h^P[ X2 s in(cjP) + y 2 cos(cjP)], (35) 
where 



G(P) = (xi - -)e^ p + 



7 



7 



+ 



a 2 + uj 1 



- 7 p 



(ay 2 — ujx 2 )(e aP cos(uP) — 1) + (uoy 2 + ax 2 )e ay sin(wP3;f|5) 



p , 



The time from one pulse to the next is determined when y 1 — with t 2 
Thus, from Eq. (28) we have a condition to determine P as 



j G(t)dt 



0. (37) 
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Finally, for the inversion of laser- 1 we have 
xi(t ) >-> [xi(ti) = G(P)] ^ xi(t 2 ), 



(38) 



yielding 

Xl ^ -G{P) + \tbG(P) 2 . (39) 

The map is evaluated as follows: (i) The current state of the system, given 
by xi, X2 and y 2 , is known, (ii) Compute the time P of the next pulse using 
Eq. (37). (iii) With P fixed we can evaluate G(P) in Eq. (36). (iv) The cur- 
rent state of the system and G(P) determine new values for x 2 and y 2 with 
Eqs. (35). (v) Finally, X\ is found from Eq. (39). Summarizing, we have 



(ii) Eq. (37) f 1 (P;xi,x 2 ,y 2 ) = 0, 

(iv) Eq. (35) x 2 ^ f 2 (x 2 ,y 2 ,G(P)) 

V2 i-> h{x 2 ,y 2 ,G{P)) 

(v) Eq. (39) xi » h{G{P)). 



(40) 



4-4 Periodic solutions as fixed points 



Fixed points of the map described by Eqs. (40) correspond to periodic solutions 
of the original flow, Eqs. (6). However, it is not feasible to analyze the map 
without further approximations. We will look for fixed points making use of 
e< 1. With heavy use of symbolic computation, we find that the maximum 
amplitudes and the period of the oscillations are given by 



max[xi] = 7T + y|^iN , 



maxfa^] = max[y 2 ] = W^ 2 — — 



3a 2 1^2 1 ' 

P = 2max[xi]. (41) 

For each result the neglected terms are 0(e). In addition, we also obtain the 
phase relationship result that x 2 ~ when Xi is at its minimum, which is 
consistent with Fig. A. 2. We have plotted the predictions of Eqs. (41) along 
with the results from numerical simulations in Fig. A. 8 and they show good 
agreement, where for clarity we have removed the higher bifurcation branches 
present in Fig. A.l. 
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The period P and the amplitude of laser-2 show excellent agreement. We see 
that max [£2] ~ l/l^l 1 ^ 2 - This may initially seem counter intuitive because the 
pulses of laser-1 grow with \8 2 \ and provide a greater kick to laser-2. Indeed, 
a leading-order approximation to the kick applied by laser-1 5i2G(P) ~ SiP, 
thus, the strength of the kick increases as the period increases. However, with 
longer periods the exponential decay due to damping in the outer regime has 
more time to decrease the amplitude of laser-2. The net effect is a decrease in 
max[i 2 ] with increasing \5i\. 

The net coupling strength of laser-2 to laser-1, with respect to \5 2 \, is 5 2 y 2 = 

1 /2 

0{5 2 )• Thus, the amplitude of the pulsations increases with increasing S 2 . 
The fit between the analysis and numerics is not as good for laser-1. To achieve 
a better fit we need to derive a map that includes higher-order terms in e, which 
we have not attempted. 



5 Discussion 

For two coupled lasers we have studied the bifurcations that occur when the 
coupling constants are of opposite sign and unequal. Specifically, the coupling 
is asymmetric in that we fix one coupling constant (5\ > 0) to be small, 
while varying the other (5 2 < 0). There are two Hopf bifurcations to periodic 
output, one for 5 2 positive and one for 5 2 negative. We have focused our 
attention on the latter because of its similarity to our work with delay coupling 
in [1]. When the output of laser-2 is nearly harmonic, the negative coupling 
effectively corresponds to phase shift by half of a period. This is equivalent to 
delayed coupling when the delay is half the period. 

As \S 2 1 is increased there is an initial Hopf bifurcation from the laser's CW 
steady-state to periodic solutions. We then observe two resonance regimes 
where the coupled system shows novel and interesting output, (i) Close to 
the Hopf bifurcation a resonance can occur where the amplitude of the laser 
oscillations becomes large. This is unexpected because both coupling constants 
are still small. The resonance is due to the negative-coupling that effectively 
reduces the damping in the laser, (ii) As the strength of the coupling increases 
further the periodic solutions may, depending on the parameters, exhibit a 
period-doubling sequence to chaos as well as the coexistence of subharmonic 
solutions. These effects are reminiscent of a periodically modulated laser. In 
the case of the coupled lasers the laser receiving the weak coupling remains a 
nearly harmonic oscillator that excites the strong resonances of the pulsating 
laser. 

The large-amplitude resonance that occurs for small coupling can be eas- 
ily understood from the well-known coupled Stuart-Landau equations given 
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generically by Eq. (21). Steady-state solutions of Eq. (21) correspond to the 
amplitude of the periodic solutions. Simple algebra shows that tuning some 
of the coupling parameters to be of opposite sign can lead to a vanishing 
denominator. Physically, the coupling term is providing an effective negative- 
damping that cancels with the lasers self-damping and, hence, a resonance 
effect. Our analysis assumed that both coupling constants were of the same 
relative size, Sj = 0(e). However, other scalings satisfy the Hopf condition, 
e.g., 5i = 0(e 1 / 2 ) and 5 2 = 0(e 3 / 2 . This does not change the qualitative prop- 
erties of the bifurcating periodic solutions in any way. 

In App. A we have attempted to describe the solutions that occur near the 
peak of the resonance. In this regime, both lasers show approximately equal 
amplitude pulsating solutions as exhibited by Fig. A.4(al). Our analysis re- 
produces the equation for the location of the resonance 62s- It also predicts 
that the period should be twice the maximum amplitude of the inversion. That 
both the period and the amplitude show a resonance peak at 62s in Fig. A. 4a 
& b is consistent with this result but the scale factor of 2 is not correct. Also, 
we do not obtain and expression for how the period (or amplitude) depends 
on the parameter 82- A difficult higher order analysis would be required to 
remedy these last two limitations. 

In Fig. A. 2 we showed that when \S 2 \ = 0(1) (see Fig. A.l) that one of the 
lasers amplitudes is large while the other's remains small. The term "local- 
ized solutions" has been used to describe the case when identical oscillators 
in a coupled system exhibit amplitudes of different scales. In coupled lasers 
localized solutions have been described by Kuske and Erneux [15] who de- 
rived a similar pair of integral conditions to Eq. (A.l). Instead of looking 
for pulsating solutions, they considered 0(1) solutions approximated using 
a Poincare-Lindstedt method, and small amplitude solutions approximated 
with the method of multiple scales. Repeating this analysis for our problem 
reproduces the Hopf bifurcation results that we obtained in Sec. 3. 

To describe the system's output when the coupling is strong, we have derived 
a map that predicts the period, amplitude and phase of the lasers from one 
pulse of laser- 1 to the next. Constructing the map relies on combining both 
strongly and weakly nonlinear asymptotic methods. That is, we used matched 
asymptotics to describe the pulsating laser- 1 and to separate one period into 
an inner and outer subintervals. For the small-amplitude laser-2 we used the 
multiple scale methods within each subinterval. We obtain very good agree- 
ment between the amplitude of laser-2 and the overall period. The amplitude 
of the pulsations of laser-1 are not described quite as well. This could be 
because we need to consider higher-order terms in our solutions, or we are 
comparing our numerical and analytical results in a less than ideal parameter 
regime. 
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The large-amplitude solutions in the resonance peak just after the Hopf bi- 
furcation are not the same as those that appear due to a "singular-Hopf bi- 
furcation" [17,16]. In the latter case, the large amplitude oscillations are due 
to crossing a separatrix separating small-amplitude solutions near the Hopf 
bifurcation from large-amplitude relaxation oscillations formed around a slow 
manifold. The functional form of the dissipation terms in the present problem 
disallows this type of behavior. 

To our knowledge, the small-coupling resonance peak has not been previously 
described; most likely this due to consideration of physical systems where 
controlling the sign of the coupling is not possible. However, in a forthcom- 
ing study we will show that for same-sign coupling but with delay, we can 
again produce the resonance because the delay provides the phase shift that 
effectively leads to the sign change. 



A Pulsating solutions for weak coupling 

In Sec. 3 we looked for small-amplitude solutions near the Hopf bifurcation 
point. We now allow the amplitude of the solution to be arbitrary but will still 
consider the coupling to be small. The laser system Eq. (6) can be rewritten so 
that the intensity and inversion evolve according to a perturbed-Hamiltonian 
system [14]. From the coupled-Hamiltonian systems, we derive solvability con- 
ditions for T-periodic solutions as 



The integrals in Eq. (A.l) are computed by evaluating Xj and yj on periodic 
orbits of the e = 0, Hamiltonian system. Unfortunately, we do not have closed 
form analytical solutions for Xj and i/j. However, for pulsating output we can 
construct approximate solutions to the Hamiltonian system using matched 
asymptotic expansions similar to what we did in Sec. (4). In this case, we 
match the outer and inter solutions to determine a uniform solution that can 
be used to evaluate the integrals in Eq. (A.l). Because we have carried out 
similar calculations in the past [4,14], we will only summarize the details of 
the intermediate steps. 

Before proceeding, we mention that Kuske and Erneux [15] derived an almost 
equivalent pair of solvability conditions for two coupled lasers. Their goal 
was to investigate so-called "localized" solutions where one laser has 0(1) 
amplitude oscillations, while the other has small oscillations, and both are 
approximated using the Poincare'-Linstedt perturbation method. Doing this 



T 









17 



calculation for our problem effectively reproduces our earlier results obtained 
near the Hopf bifurcation point and thus does not provide new information. 

We assume that laser- j has period Tj, which is some fraction of the total 
period with T = nfTj. The first term in each integral is x 2 - and because it does 
not involve the other laser is independent of the phase relationship between 
lasers j and k. Thus, using the results from [14] we have 

]x]dt = n fTl (A.2) 
o 



Integrating Xji/k is more complicated because we must allow for a phase differ- 
ence, 7^, between the two pulsating lasers. However, it is easy to predict the 
form of the result. The intensity yt is pulsating and acts like a delta function 
that samples the inversion Xj at the time of the pulse. The effect of the inte- 
gral is to sum all of the sample values of the population inversion. In effect, 
we have a pulse train due to one laser sampling the population inversion of 
the other. After carrying out the detailed calculations based on the approxi- 
mate solutions of the Hamiltonian system, we obtain our final result for both 
solvability conditions: 

- a^Tf + T 2 5 2 g + 2)T 2 ) = 0, (A.3) 



where 



-jTi + Tf + t, m-T^ <t<(j + |)7\ — Ttj, 



-(j + 1)71 + T 4> + t, (j + -T 4> <t<(j + 1)71 - T 



and 



n n^-l -I 

a'Tl + Ti5\ x 2 ((j + -)T, - T+) = 0, (A.5) 
1Z j=0 z 



where 



, -kT 2 + t, kT 2 <t<{k + \)T 2 

x 2 (t) = { . (A.6) 

-(k + l)T 2 + t,(k + \)T 2 <t<(k + l)T 2 



The inversion variable Xj of each laser is a saw-toothed type function that 
increases linearly from the time of the previous pulse to the next. Specifically, 
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for laser-2, x 2 increases from (at time kT 2 ) to x 2 = T 2 /2. The intensity pulse 
depletes the inversion to x 2 = — T 2 /2, whereupon x 2 then increases linearly 
back to 0. Laser-1 is the same except that we must allow for a phase time 
between the two lasers. 

We consider the simple case of a 1:1 resonance between the lasers where T = 
T\ = T 2 so that n\ — n 2 — 1. The solvability conditions reduce to 

-^T 2 + fe(|) = (A.7) 
-^T 2 + te(|-T,) = (A.8) 



Because we are interested in periodic solutions, it is sufficient to consider 
< T,!, < T. Then, substituting for x\ and x 2 , the conditions reduce to 

-^T 2 + 5 2 (-| + T,) = (A.9) 
-^ + 5^-^ = (A.10) 



After eliminating the phase T^, we obtain 

(1 + ^)T 2 = 0. (A.ll) 

CL\0\ 



For periodic solutions with T ^ 0, we are forced to set the term in parenthesis 
equal to zero. This is exactly the same condition that identifies the location 
of the singularity in the Hopf bifurcation equation (20). This confirms that 
there is an equal-amplitude, pulsating 1:1 resonance between the lasers when 
$2 — $2S- However, we do not have any information on the period or amplitude, 
which would require continuing the analysis to higher order. 
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Fig. A.l. Numerical bifurcation diagram (AUTO) [18] using 82 as the bifurcation 
parameter. Note, the coupling is negative so 82 < 0; as we increase l^li the coupling 
is increasingly negative. Thick (thin) lines indicate that the periodic solutions are 
stable (unstable). There is an initial Hopf bifurcation (HB) from the steady state 
to periodic solutions. As l^l is increased the initial periodic orbit exhibits a pe- 
riod-doubling (PD) sequence of bifurcations to chaos; we show only the first two 
PD bifurcations. Simultaneously, subharmonic periodic solutions appear through 
saddle-node (SN) bifurcations and will also period double. For these parameter val- 
ues, {82s I > \$2h\ so that the small-coupling feedback resonance peak is exhibited. 
(Fixed parameters are e = 0.001, b\ = 62 = 1, Q>i = «2 = 25 and 8\ = 0.04.) 
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Fig. A. 2. For l^l = 2.18 (62 < 0). (a) Laser-1 has pulsating intensity (solid) and a 
triangular-shaped population inversion (dashed) because it is strongly modulated 
by laser-2. (b) Laser-2 is nearly harmonic because it receives only weak coupling 
from laser-1. In (a), the times marked to, t\ and £2 define the outer and inner regimes 
discussed in Sec. 4. 
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Fig. A. 3. Comparison of the amplitude of the population inversion of laser-1 (left 
axis, solid curve) and the inversion of laser-2 (right axis, dashed line) for: (a) Before 
the first period-doubling bifurcation, l^l = 2.18 (same as Fig. A. 2). (b) After the 
period-doubling bifurcation 1^1 = 5.04. (c) Just after the saddle-node bifurcation 
near the limit point, j^l = 0.70. There are three maximums of laser-2 for every two 
maximum of laser-1. 
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Fig. A. 4. (a) Amplitude of inversion X2 as a function of 82 = ed,2 after the Hopf 
bifurcation; numerical (+) (Auto [18]), analytical from Eq. (17) (solid line). Param- 
eter values are e = 0.001, a\ = a2 = 2, d\ = 3 and a = so that d,2H = 4/3 and 
(I2S = 3. In the inset (al) we show the pulsating intensity near the peak of the 
resonance (yi dashed. 1/2 solid), while in inset (a2), the intensity has returned to 
be small-amplitude and nearly harmonic, (b) Period of oscillations as a function of 
62 = ed.2 with analytical result (solid line) from Eq. 19. 
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Fig. A. 5. Amplitude of X2 as a function of 82 = ed.2 after the Hopf bifurcation; 
numerical (+) [18], analytical (solid line). The parameters have been chosen as 
02 = 2.9 and d% = 3.1 so that the singular point d,2s = 2.13 is very near the Hopf 
bifurcation d,2u = 1.87. 
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Fig. A. 6. Amplitudes for the case of three identical lasers (fixed parameters are 
cti = 02 = «3 = 2, b = 1, Pi = fa = 1) e = 0.001.) The fixed-coupling contants 
are 613 = 623 = 1.3 and 621 = #31 =3. The coupling of laser-2 into laser-3 is 
positive with size £32 = l^l, while the coupling of laser-2 into laser 1 is negative 
with 612 = 62 < 0. 
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Fig. A. 7. Same data as Fig. A. 6 but with the amplitudes of laser-1 and laser-3 shown 
simultaneously. The solid curve shows the stable solutions that appear after the Hopf 
bifurcation (HB). There is a stable resonance peak for \d,2\ ~ 3. The dashed curve 
shows unstable solutions that exist between the two saddle-node (SN) bifurcations. 
There is an unstable resonance peak for \d,2\ ~ 1.75. Both laser-3 and laser-1 (and 
laser 2) exhibit large amplitudes during the unstable resonance, while only laser-1 
(and laser-2) has large amplitude during the stable resonance. The dotted curve is a 
projection of the actual data into only the laser-1 plane to compare to laser-1 data 
in Fig. A. 6. There is bi-stability between the saddle-node (SN) points. 
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Fig. A. 8. Comparison of numerical bifurcation results (solid) to the analytical results 
(dashed) from the map in Sec. 4, specifically, Eqs. (41). Parameters are the same 
as in Fig. A.l. 
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